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Abstract 

Two different types of transitions of the superfluid ^He to quantum turbulence regimes are 
studied for ID geometry in the case when the influence of the normal fluid on superfluid flow is 
suppressed. It is shown that the roton mechanism of transition to quantum turbulence leads to 
a critical velocity satisfying the relation Vc oc d~^^^. In the super-flow mechanism, the transition 
to quantum turbulence arises when the "quantum Reynolds number" is about 10'^ and the critical 
velocity depends on channel size d as Vc cc d^^ in agreement with the equations of motion for a 
superfluid component of the liquid ^He being disturbed by small fluctuations of the normal fluid. 

PACS numbers: 67.25.dg, 67.25.dj, 67.25.dk 



1 



The understanding of the existence of critical velocity in superfluid ^He and its micro- 
scopic nature is a long-standing problem which also closely related to superfluid turbulence 

O]. The critical velocity in a superfluid is the threshold above which the flow of the 
superfluid component becomes dissipative and the property of superfluidity is lost. The first 
well known criterion for critical velocity, Vc = min e{p)/p where e{p) is the energy spectrum 
of the elementary excitations at momentum p, was proposed by Landau This criterion 
yields the value Vc — QO m/ s which does not depend on capillary parameters. The Feynman 
criterion ^, following Onsager, is based on the quantization of superfluid circulation and 
defines the velocity Vc — {h/md)ln{d/a) at which the vortices can be excited. Here d is the 
channel size and a ~ 4A is the vortex core radius. These two different criteria can differ 
greatly from experimental data [5]; however, the Feynman estimation is much closer to ex- 
perimental values. The difficulty in this problem also arises due to the existence of different 
breakdown mechanisms of the superfluidity 5|, ISl, lZJ . Among the several problems associated 
with critical velocity the type of elementary excitation which would be responsible for the 
onset of dissipation and for a critical velocity in the superfluid plays an important role. In 
Feynman's approach such elementary excitations are vortices excited in the superflow of the 
channel [J]. 

We study in this letter the existence of the critical velocity in superfluid ^He when the 
normal fluid does not influence the superfluid flow. The most efficient method of preventing 
the normal fluid from interfering with the flow is by introducing superleaks at the ends of 
the capillary js]. The narrow paths in the fine powder act as solid walls to the normal fluid 
and prevent it from entering the flow region. As soon as the normal component is sufficiently 
suppressed by application of superleaks in the He II flow path between parallel plates or in 
a capillary it is found that the experimental results are well described by the temperature 
independent critical velocity Vc oc c?"^/^ jsl, [lol, 11, I21, 13] when T < 2 K and the separation 



of the plates is do < d < dc where do = 0.8 ■ 10 ^ cm and dc = 1 cm. 



Our explanation of this relation for Vc is based on the cluster model of the roton 

mm 

in liquid ^He. It was shown in jl^-fisjl that the roton is an excitation of a cluster of 13 atoms 
having a central atom surrounded by a shell of 12 atoms situated at the vertices of a regular 
icosahedron, although in this picture each atom can deviate about its location in the cluster. 
We show that rotational excitations of the roton cluster of 13 atoms in the superflow are 
the elementary excitations responsible for the onset of dissipation and hence for a critical 
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velocity in the superfluid for a range of capillary widths d. 

The elementary rotational excitations are excited by collision with the surface of the 
channel. The rotons can exchange between each other by the rotational energy in the 
superfluid and when the velocity of the flow becomes greater of some critical value Vc the 
number of the elementary rotational excitations grows in the form of an avalanche and the 
superfluid motion breaks down because the irrotational superfluidity requirement curlvg = 
for macroscopic superfluid velocity is no longer satisfied. In this microscopic picture we 
assume that interaction of these rotational excitations forms the vortices in the superflow 
of the channel. It is shown below by dimensional analysis that the roton mechanism of 
transition to quantum turbulence leads to a critical velocity given by Vc oc d~^^^, in agreement 



with experimental data 



0, Q, Q, 



12, 



13|. 



We also consider in the paper another type of turbulence (super-flow quantum turbulence) 
which is quasi-classical and arises when the size between parallel plates d > 1 cm and the 
dimensionless parameter TZ = d/X is about 10'^ where A = 2TTh/mv is the characteristic 
de Broglie wavelength of helium atoms in the channel flow. It is found that a super-flow 
quantum turbulence arises when v > Vc where the critical velocity satisfles the relation 
Vc OC d^^. This type of turbulence differs considerably from the quantum turbulence when 
the normal fluid is not suppressed [l, 2] and, in particular, the critical velocity by Feynman 
criterion in this case leads to a value about three orders smaller than was found in the 
experiments jisl at (i ^ 1 cm for normal mass density p„ ^ p. 

The Schrodinger-type equation for the wavefunction \I'(ai, 02, as, t) describing the roton 
cluster in the space of ellipsoidal parameters {k = 1, 2, 3) as a complete system of Nc 
bound atoms is Q, Q : 

= [-2M, £ ^ + ^ £ "^^^ + ^ j ^' 

where Mq is the effective mass of the roton cluster and Q is an interaction constant given by 

Mo = -^, G = . (2) 

For spherically symmetric roton states (a/c = a) the frequency uj = uj^ (fc = 1, 2, 3) is 



h /l5(iV,-l)ao 
^ = — = ' (3) 
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where a is the effective radius of the roton cluster and = 2.2 A is the s-scattering wave- 
length of ^He atoms: 

ao = — U{r)r^dr. (4) 



V 47rp y h' 

Here p is the mass density of liquid helium and the s-scattering wavelength is calculated by 
intermolecular Lennard- Jones potential U{r) = 4e[(ro/r)^^ — (ro/r)^] for ''He. The stability 



condition of the roton cluster has the form Ad > 2a where Xd = 2-Kfi/\^3mkBT is the 
thermal wavelength. Using the quadratic approximation in the potential of Eq. ([T]) the 



eigenenergies 



14 



15| to good accuracy are 
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En,n2m = Vo + J2^^k{nk + l/2). (5) 

k=l 

Here rik = 0, 1, 2, ... (at k = 1, 2, 3) are the quantum numbers of a 3D quantum harmonic 
oscillator. In the spherically symmetric case Eq. ([1]) yields the eigenfrequencies Qk as 

n^ = Q^ = y/2uj, = V5lj. (6) 

In the general case rotational quantum states of the roton clusters can also be excited. If an 
inequality ^ = Bc/huo ^ 1 is satisfied where = Ti' jll^ is the rotational constant and Jc is 
the moment of inertia of the roton cluster, then for small vibrational quantum numbers 
and small angular momentum quantum numbers J the eigenenergies of the roton are 

Kn,n, = ^0 + E ri^kiUk + 1/2) + BJ{J + 1). (7) 
k=\ 

This equation comes from standard quantum mechanical methods in Eq.([T]) as developed 
for rotational-vibrational molecular spectra. In this general case the number of atoms in 
the roton cluster satisfies the condition [l^ 

^i^n^no, - -^000 = ^> A = e - eo, (8) 

where A is the roton gap defined by the parameters e = —U{rra) and eo = — f/(2ao). Here 
Tm = 2'/^ro is the distance for which the intermolecular potential energy U (r) has a minimum 
and 2ao is the mean distance between helium atoms in the liquid *^He. A self-consistent-field 



Hartree-Fock method 



16|, 



ITj l yields the parameters: e/ks = 10.6 K and = 2.98 A, hence 
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from Eq.([8]) it follows that t^jhs = 8.65 K. Using Eqs. ([MH]) the equation for the number 
Nc in the roton cluster is 



V2{ni + na) + VSrigl + BJ{J + 1) = A. (9) 



Because we assume ^ <^ 1, the term BcJ{J + 1) can be neglected and from Eqs. 
follows that the integer number A^^ is given by 

- a[V2{n, + ns) + V^n^fiN, - if = 0, (10) 

where a = (Aq/A)^, and the parameter Aq is 



^0 



^o=M- ■ (11) 



The smallest solution to Eq. (|TOl) is Nc = 13.1 fisl] in a symmetric vibrational state with 
Hi = n2 = and n3 = 1, hence the most stable roton clusters should be those consisting of 13 
helium atoms. We note that for roton clusters with A^^^ = 13, ^ = (9(10^^) and a = 5.22 A. 

The kinetic energy flow density of the superfiuid between parallel plates for the case when 
there is no normal fluid in superfiuid flow (p„ = 0) is 

„ 1 1 O 1 q „\ 

where S = Ld is the area of the cross section and p = ps- Let us consider ID geometry (the 
width of the plates L is much greater then the distance d between plates); then the kinetic 
energy flow per unit length along the width of the parallel plates is given by / = Fd or in 
the limit L oo: 

The function / has such a form only for Vg < fc because when the velocity of the flow is 
greater then Vc the onset of dissipation occurs in the liquid ^He. Hence the critical value of 
/ given by fc = pv^d/2 can be treated as a threshold parameter of the superfluid. Another 
important parameter of the superfluid is the minimal rotational energy Sr = Eqq-^^—Eqq^ which 
rotons can exchange in the superfluid between the rotational states with angular momentum 
quantum numbers J = 1 and J = 0. From Eq. (J?!) it follows that this parameter is deflned 
by the moment of inertia of the roton cluster Ic as Sr = 2Bc = fi^/Ic- 

We may assume that when the normal fluid is suppressed the critical velocity depends 
on three parameters: the threshold parameter fc, the minimal rotational energy Sr and the 



minimal angular momentum h which rotons can exchange in the superfluid. The parameters 
Er and h are connected with the dissipation process of the energy and angular momentum 
in the liquid "^He at Vg = Vc- Thus it is assumed that Vc = Q{fc,£r,^) where Q is some 
unknown function of three variables. Taking into account that the threshold parameter has 
the form = pv^d/2 the relation for critical velocity can be rewritten in an equivalent form 

Vc= Q{pd,er,h), (14) 

where Q is a new unknown function of three independent variables pd, Sr and h. It can be 
shown by dimensional analysis that there exists a unique function Q of the variables pd, 
and h which has the dimensionality cm s~^. This function has the form Q = '-){pd)^e^h^ 
where 7 is a dimensionless parameter and the unique powers are: n = —1/4, k = 3/4 and 
/ = —1/2. Hence, the general relation given by Eq. (fT^ yields the unique equation 

Vc = rd-'/\ (15) 

where the constant parameter F has an explicit form: 



We emphasize that the Eq. ( ITSl) coincides with the known empirical relation [9|,lld,lll|,ll2|,ll3| 
for critical velocity in the case when the normal component is sufficiently suppressed by 
application of superleaks. Moreover, the explicit expression given by the Eq. ffT6|) a llows 
us to estimate the constant parameter T. It is assumed in dimensional analysis [iSj that 
the dimensionless parameter (in our case 7) has a magnitude of order unity. More exactly, 
one may guess that (9(10^^) < 7 < (9(10). The evaluation of the constant parameter F 



by Eq. ( |T6l) for roton clusters with Nc = 13 and 7 = (9(10 ^) yields T = O 
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}cm 



5/4.-1 
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13| 



which is in agreement with the experimental value F = 1 cm^/^s~^ in Ref. 
Because there is no other characteristic energy in liquid ^He comparable to Er in size, this 
estimation demonstrates that roton breakdown mechanism of the superfluidity explains the 
experimental relation Vc = Td~^^'^ for critical velocities when do < d < d^- We also emphasize 
that the Eqs. f[T^ 16) shows that the superfluidity is a quantum phenomenon because Vc —>■ 
in the classical limit h —>■ 0. 

The experiment also demonstrated the temperature dependence of the critical velocity 
[isl showing that Vc dropped suddenly to zero above about temperature Tq = 2 K and below 



this temperature Vc was constant. We may estimate the characteristic temperature Tq by 
condition Xd > 2(i// where dn = 2ao = 4.4 A is the effective diameter of hehum atoms. This 
condition means that the superffuidity takes place only for temperatures T < Tq and leads to 
characteristic temperature: Tq = -k'^U^ /VlmkBO^- This yields the temperature Tq = 2.05 K 
which coincides with the experimental characteristic temperature. 

We assume that for thicknesses of the channel d > dc where dc is the critical value for 
existing of the relation Vc oc d~^^^, the physical breakdown mechanism of the superfluidity is 
different and it is connected with the characteristic interaction lengths I > 2d = 10.44 A. In 
this case the wavenumbers belongs in the phonon region k < ko where ko = n/d = 0.6 A^^. 
It can be shown that in this region of wavenumbers the effective pair interacting potential 
has the form (which will be published elsewhere): f/(r) = U(r) for re [aQ, a^] and U{r) = 



for r ^ [accim] where U{r) is the Lennard- Jones potential [l6|, llTI] and = 2rm — ao = 
3.76 A. This means that in the effective Hamiltonian for liquid ^He there is a cut-off in the 
intermolecular potential for distances r < and r > a^. The first cut-off r < ao is related 
with the repulsive part of the intermolecular potential in the modified Born approximation 
[is! ] (see Eq. (jl])). The cut-off in the region r > am, where am is defined by relation 
fm — (^0 = o^m — fm, IS couuectcd with the long-range attractive forces between one atom and 
atoms in the surrounding bulk liquid summing to zero. 

The many-body effective Hamiltonian for the region of wavenumbers k < k^ yields the 
Heisenberg equation of motion for the field operator as 

at 2m 

where the coupling constant G = Arc U{r)r'^dr = 1.79 ■ 10"^"^ g cm^ s~^. We note that 
the parameter G significantly differs from the coupling parameter g = Anaoh'^ /m = 4.63 • 
-|^g-38 g g-2 Gross-Pitaevskii equation for dilute helium gas [2^. Using the field 
operator in the form 'ip{r,t) = ipsi^^t) + '?/;„(r, t) where {■ip{Y,t)) = ips{v,t) and {i^n{'r,t)) = 
we may write the full mass density of the liquid helium as the sum p = ps + pn with 

ps{r,t) = m|V's(r,^)|^ pn{r,t) = m(V'l(r, t)V^„(r, t)). 

We consider the case pn Ps', then in the dynamic Popov approximation [2]J the anomalous 
density {i^ni^n) and the three-field correlation function {'ii}\'tpn'>Pn) can be neglected in an 
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FIG. 1: Flow velocity v for ^He in a channel of width d showing regions of superfluidity and 
turbulence. This is for ID geometry with T < 2K and the influence of the normal fluid on 
superfluid flow suppressed by superleaks. The line shows the critical velocity Vc which separates 
superfluidity from quantum rotational turbulence for widths d < 1cm (Eqs. (115^ 16)) and from 
quantum super-flow turbuluence for widths d > 1 cm (Eqs. ()221 23)). 



equation for macroscopic wavefunction ilJs{r,t): 



ih 



dips 



h' 



(17) 



dt 2m 

The Heisenberg equation also yields the linearized equation for the field operator ^^^(r, t) as 



ih 



dt 



h^ . . 

2m 



(18) 

tpsoGxp^—ifit/h) where 
^so = \/Pso/'"^exp(i$o) is the constant 
amplitude. In this case the solution of Eq. ( fT8l) for the field operator '?/'„(r, t) is 

JT^w ' 



The homogeneous solution of Eq. f|T7|) at p„ ^ ps is 
fi = Glipsol"^ is the effective chemical potential and ips 



where fx = mc^, 7fc = e '^^^°{huj{k) — h^k!^ /2m — mc^)/mc^, Ck, are the annihilation and 
creation Bose operators for excitations in hquid ^He and uj{k) is given by: 



-(^) = ^f^ + 5^- (19) 

Here c = ^/Gp^/m is the velocity of sound in He II, p^o = p[l ~ (T /Tx)^-^] is the density of 
the superfluid component of the liquid ^He assuming thermodynamic equilibrium 191], and 



Tx = 2.18 K is the critical temperature. As an example for T = 1.1 K, c = 2.4 ■ 10^ cm/s. 
The Bogoliubov's type spectrum in Eq. flTI?]) for acoustic waves excited around the constant 
density pso at wavenumbers k < k^is close to the linear spectrum uj{k) = ck coinciding with 



the experimental data for the phonon spectrum in He II [22(] when T <^T\ {pn ^ Ps)- This 
solution for the field operator ■?/'„(r,t) yields the effective Hamiltonian in the diagonal form 
H = Eq + J2k fi'U!{k)cl^cii which at wavenumbers |k| < k^ describes the phonon excitations. 
The Eq. ( |T7I) can be written in typical hydrodynamic form in terms of the superfluid 



density ps and velocity Vg by standard definitions ips = \J Ps/mex.\){i^) and Vg = 

(^/m)V$. Using the characteristic size d between parallel plates (ID geometry) and the 
flow velocity f , new dimensionless variables can be defined by t' = vt / d, r' = r / d, = Vg/f, 
p'g = Gps/m'^v'^, p^ = Gpn/Tn^v"^ and the dimensionless hydrodynamic form of Eq. (ITTI) is 

+ (v^ ■ V')v^ + V'pI = i^V'J^'] - 2Vp'^, (20) 



dt' ^ ' ' ; s . ' g^27^2 

^ + ^\p's<) = 0. (21) 

Here TZ = mvd/2Tih is the dimensionless number characterizing the flow of the superfluid 
^He in the channel. Using Eq. 0181) one obtains further equations for p'„. These are found 
also to depend only on IZ and no other parameters. Thus IZ, the "quantum Reynolds 
number" is the unique dimensionless variable parameter completely characterizing the flow. 
We note that if both terms on the right hand side of Eq. (120|) are negligible (Thomas- Fermi 
approximation) then we get the classical Euler equation for potential flow of a nonviscous 
fluid with pressure P = Gp1/2m?. For large numbers TZ the quantum pressure term [l^ 
which scales as TZ~^ in Eq. (120]) becomes less important in comparison with the non-linear 
inertial term (Vg ■ V')Vg and the normal fluid fluctuation term — 2V'p[j. In this case, as it 
is observed experimentally, for some large 7^ (at c? = 1 cm and v = 1 cm s~^) the laminar 
flow becomes increasingly unstable leading eventually to turbulence [l3| . Because the type 



of flow depends only on the unique dimensionless number TZ one may conclude that the 
quantum super-flow turbulence (QSFT) arises when TZ > TZc where TZc is critical number. 
This condition yields the QSFT sX v > where the critical velocity is 

2Txhnr 



md 



(22) 



The critical distance dc is the boundary point on the curves given by Eq. (fT5i) and Eq. (1221) . 
which yields the critical number TZc as 



mVdl'^ 



The experimental value 



13| of the critical width is = 1 cm and hence from Eq. (1231) it 



follows that TZc — 1000. This large number means that the de Broglie wavelength A of the 
helium atoms in the channel flow is much less then the channel size d for TZ = TZc- Thus 
the QSFT is the quasi-classical regime because the limit TZ ^ oo formally is equivalent to 
the limit % 0. However we note that this regime arises from fluctuations of the normal 
component (see Eq. ( fTSi) ) and has a quantum nature because TZ is finite. Fig. 1 shows 
the transition from superfluidity to QRT and QSFT regimes as predicted by Eq. ( JT5l) for 
d < 1cm [vc = d~^/^ c.g.s.) and Eq. ( l22i) for d > 1cm (vc = d~^ c.g.s.) respectively. 

In conclusion, in this paper are considered two different mechanisms of transitions of 
the superfluid ^He for ID geometry to quantum turbulence regimes when the influence 
of the normal fluid on superfluid flow is suppressed by superleaks. The first mechanism 
takes place for thicknesses 0.8 ■ 10"^ cm < < 1 cm and T < 2 K, and can be explained 
by excitation of the rotational states of the rotons in the vicinity of the surface of the 
channel and further exchange of the rotational energy between roton clusters. In this case 
the functional form of the critical velocity is given by Eq. (fT5|) and matches the observed 
trend in experiment. The second mechanism of transition to quantum turbulence arises 
at > 1 cm and TZ = TZc — 1000 when the quantum pressure term is diminished and 
the laminar flow becomes unstable. In this case the functional form of the critical velocity 
is given by Eq. ( l22l) — a prediction whose validity (of interest) requires the extension of 
existing experimental data. 

I am grateful to Dr M. J. Collett and particularly to Dr D. Wardle for numerous valuable 
comments and useful discussions of this work. 
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